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We report our studies of the potential energy surface (PES) of selected binary Lennard- Jones 
clusters. The effect of adding selected impurity atoms to a homogeneous cluster is explored. Inherent 
structures and transition states are found by combination of conjugate gradient and eigenvector- 
following methods while the topography of the PES is mapped with the help of a disconnectivity 
analysis. We show that we can controllably induce new structures as well as reorder and stabilize 
existing structures that are characteristic of higher-lying minima. 
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I. INTRODUCTION 

The minimization/optimization problem is one of 
the more ubiquitous and challenging in computational 
science^. Central to researchers in the physical sciences 
and engineering, this problem is also of primary impor- 
tance to social, biological, and economics investigators. 

Driven in large measure by such widespread interest, 
there has been appreciable progress on the minimiza- 
tion problem. Especially notable have been algorith- 
mic advances in the form of annealing and stochastic 
relaxation appro achesiSi^i^i^ as well as basin-hopping 
techniques''^. In both classical and quantum form, these 
methods offer valuable, complementary alternatives to 
traditional, gradient or pseudo-gradient approaches^. 

In addition to algorithmic developments relevant to the 
minimization problem, there have also been notable ad- 
vances in the tools to classify and analyze the topogra- 
phy of the underlying objective functions. In chemical 
applications, the principal focus of the remainder of our 
discussion, the objective function of interest is typically a 
specified potential or free energy surface. Following Still- 
inger and Weber^'^^, it is useful to perform an "inherent 
structure" decomposition of the associated configuration 
space by employing the minima (local and global) of this 
surface. These inherent structures, their relative order- 
ings, and their connectivity provide important informa- 
tion concerning the structure, function, and dynamics of 
the associated physical system. Disconnectivity analy- 
sis introduced by Czerminski and Elber^^, discussed by 
Becker and Karplus^^ and developed by, among others, 
Wales, Doye and Millerii^ii^ii^ has proved especially valu- 
able with respect to these latter tasks. 

As evidenced by the development of classical and quan- 
tum annealing methods, there is an important inter- 
play between minimization and the Monte Carlo sam- 



pling problem. Both applications, for example, are con- 
cerned with overcoming barriers that inhibit the inter- 
conversion or isomerization of the various inherent struc- 
tures of the problem. Consequently, developments in 
one field contain implications for developments in the 
other. Advances in rare event sampling methods, such as 
J-walking^^'^' and parallel tempering^^'^^'^^'^^ methods, 
thus contain implications for the minimization problem. 

The field of atomic and molecular clusters has been 
and continues to be an important test bed for the de- 
velopment and application of minimization and analy- 
sis methods. Utilizing the methods outlined above, re- 
searchers have produced a relatively coherent picture of 
the relationship between the nature of underlying PES 
and the physical properties of the associated systems. 
For example, from the studies by Berry et a/^2i2iS^i2^, 
the single component studies of Wales et alM^^^i^^ the 
mixed LJ cluster studies of Jordan et al?^ as well as the 
research of others-^L^MMS., we have begun to understand 
the nature of systems for which the lowest inherent struc- 
tures can or cannot be readily located. 

In present paper we would like to build upon advances 
in the minimization problem by effectively turning the 
logic "upside down". That is, instead of asking what 
we have to do in order to locate or sample the global 
minimum of a specified potential energy surface, we wish 
instead to ask how we might go about controllably induc- 
ing new structures as well as reordering and stabilizing 
existing structures that are characteristic of higher-lying 
local minima. Basically, we seek to utilize what we have 
learned about what it takes to avoid local minima to in- 
stead controllably produce them. 

In principle, one can envision efforts involving both 
thermodynamic and kinetic approaches. In the present 
work we shall focus principally on the thermodynamic 
issues. Furthermore, we shall limit the discussion in the 
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present work to applications involving clusters. As dis- 
cussed elsewhere^^, clusters are of appreciable technolog- 
ical importance, are valuable as prototypes for the study 
of the properties of extended systems, illuminate issues 
related to the size-dependence of selected physical prop- 
erties, and provide valuable test beds for the development 
and application of emerging computational techniques. 
This combination of formal, computational and techno- 
logical interest has produced a vast and growing cluster 
literature2ME;21. 

The remainder of the paper is organized as follows: 
In Section m we outline the computational details of the 
present study. We discuss the methods we use to deter- 
mine the inherent structures and transition states of a 
specified cluster's potential energy surface. Using these 
methods, we examine specific results for two prototype 
systems in Section Hill These particular results are de- 
signed to demonstrate "proof of principle" with respect 
to the basic objectives of the present study for selected 
systems. Finally, in Section HVl we summarize our results 
and speculate about likely future research directions. 



results as opposed to results for particular physical sys- 
tems, it is convenient to reduce the number of free param- 
eters. To do so, we shall assume in the present study that 
the "mixed" Lennard- Jones values are obtained from the 
"like" Lennard- Jones parameters via usual combination 
rules^i 

Furthermore, we note that with the mixed Lennard- Jones 
parameters specified as in Eqs.(|3|) and Q), the inherent 
structure topography of the "reduced" potential energy 
surface of the binary system (i.e. Vtot/^xx) ^ function 
of only two parameters, (cr, e), the ratios of the corre- 
sponding adatom/core length and energy parameters 

^ = (^yy/^xx (5) 



II. COMPUTATIONAL DETAILS 

The present Section describes the computational de- 
tails of our investigations involving binary clusters of the 
form XnYm- Our overall interest will be to explore the 
extent to which we can utilize the "adatoms" (i.e. the 
Y-system) to induce, reorder and stabilize selected in- 
herent structures in the "core" X-system. While one can 
easily imagine applications involving both more and more 
complex components, we feel these relatively simple, two- 
component clusters are a convenient starting point for an 
initial study of the issues we raise. 

We shall assume in what follows that the total poten- 
tial energy is composed of a pairwise sum of Lennard- 
Jones interactions. Specifically, we assume that the total 
potential energy, Vtot, for an N-particle system is given 

by 

N 

Vtot = ^Vijirij), (1) 

where the pair interaction as a function of the distance 
between particles i and j, r^j, is given by 

%-(r'.,)=46,, [(^)i^-(^)^ (2) 

In Eq. (O the constants e^j and CFij are the energy and 
length-scale parameters for the interaction of particles i 
and j. 

For a two-component system, we must specify both 
the "like" (X-X, Y-Y) as well as the "mixed" (X-Y) in- 
teractions. With an eye toward studying trends in the 



^ = ^yy/^xx' (6) 

If necessary for a discussion of a specific physical sys- 
tem, the absolute bond lengths, energies, activation en- 
ergies, etc. can be obtained from the corresponding "re- 
duced" results by a simple rescaling with the appropriate 
core-system Lennard- Jones parameters. 

The computational task in our study is thus one of ex- 
ploring and characterizing the (reduced) potential energy 
surface of our binary cluster systems as a function of the 
number of (core, adatom) particles, (n,m), and for given 
(cr, e) ratios. In typical applications the lowest N/^ inher- 
ent structures and the associated disconnectivity graphs 
are determined. For the applications reported here, N/^ 
is generally of the order of a few hundred (thousand) or 
less. Depending on the size of the cluster, inherent struc- 
tures are found either via conjugate gradient methods 
starting from randomly chosen initial configurations, or 
by more systematic surface exploration methods such as 
those outlined by Wales and co-workers^^ and by Jordan 
et ai^^. In all cases, the inherent structures that are 
located are confirmed to be stable minima via a stan- 
dard Hessian analysis. To reduce the chance we miss 
particular local or global minima, we monitor the num- 
ber of times individual inherent structures are found and 
demand that each of the N/^ inherent structures be lo- 
cated a minimum number of times (at least 10) before we 
terminate our search. Once we are satisfied we have lo- 
cated the relevant inherent structures, transitions states 
linking these stable minima are obtained using the eigen- 
vector following methods outlined by Miller and Cerjan^^i 
and further developed by Simons et a/^i^Si^, Jordan et 
al.^^ and Wales^^. Finally, with the requisite inherent 
structures and barriers in hand, we perform a disconnec- 
tivity analysis using methods outlined by Czerminski and 
Elber^^, Becker and Karplus^^ and Wales et alM^ . 
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FIG. 1: The only stable inherent structure for X5 LJ cluster. 
Its energy (in units of the LJ well depth) is -9.104. 



III. NUMERICAL RESULTS 

In the present Section, we wish to illustrate the general 
themes we introduced in Section We do so by demon- 
strating that we can accomplish three basic objectives. 
Specifically, we show that by adding selected "impurity" 
atoms to bare "core" systems, we can: 

1. induce new "core structures" 

2. reorder the energies of existing core inherent struc- 
tures, and 

3. stabilize selected inherent structures by controlling 
the activation energies that determine their isomerization 
kinetics. 

For purposes of illustration, we shall examine numer- 
ical results for a few, simple Lennard- Jones systems in- 
volving five and seven core atoms, systems well-known 
from previous studies to have one and four energetically 
distinct inherent structures, respectively. The inherent 
structures and their associated energies for these core 
systems are illustrated in Fig. ^^nd [3 

We first consider mixed clusters of the generic type 
X5Y2. Here two impurity Y-atoms are added to the par- 
ent, five-atom X-core. We have chosen this system be- 
cause it builds upon the very simple five-atom core, a 
system that has only a single inherent structure, and be- 
cause the total system has a total of seven atoms, a magic 
number for icosahedral growth in homogeneous systems. 
Using the techniques of Section II, we then determine the 
lowest several inherent structures for a range of (cr,e) [c.f. 
Eq.© and Eq.Q]. As can be seen from Fig.|31 the total 
potential energy [Eq.Q] of the lowest inherent structure 
for the X5Y2 system shows no appreciable structure as a 
function of the (cr,e) parameters. 

On the other hand, we see in Fig. E| that the core po- 
tential energy, defined as the potential energy of inter- 
action for only the core X-atoms, of the minimum (to- 
tal) energy cluster clearly breaks into extended regions, 
each corresponding to a well-defined core structure. The 
reader should notice that each region in Fig. 21 contains 
the same "kind" of core structure but their core energies 
are slightly different. We have chosen a single "average" 
core energy value to represent all energies in the corre- 
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FIG. 2: The four, energetically distinct, stable inherent struc- 
tures for X7 LJ cluster. The energies (in units of the LJ well 
depth) are: (a) -16.505, (b) -15.935, (c) -15.593, (d) -15.533. 
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FIG. 3: Etot(cr,e) (c.f.Eqs[I] Eland EJ for the X5Y2 system. 
Note the relative lack of structure in the (cr,€) variation of the 
total cluster energy. 



sponding domain for plotting convenience. 

The distinct core structures, shown in Fig.^ have been 
identified by examining their core energies (E^ore) 
their principal moments of inertia. For each structure 
a triplet of values (Ecore, I2, I3) has been associated, 
where I2 and I3 are the moments of inertia about the 
principal axes 2 and 3, respectively. We have defined I2 
and I3 in the following way: I2=l2/Ii, 13=13/11 where I^, 
I2 and I3 are the principal moments of inertia obtained 
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by diagonalizing the inertia tensor of the system. If the 
triplet of values has not been sufficient to identify a core 
structure then we have examined the structure visually. 




FIG. 4: Ecore(cr, e) for the X5Y2 system. Here the "core" en- 
ergy is defined as that portion of the potential energy arising 
from only the core-core atom interactions. Unlike the total 
energy, the (cr, e) variation of the core cluster energy exhibits 
relatively well-defined regions. The labels of each of these re- 
gions in the figure correspond to the distinct core structures 
shown in Fig.|5] 

Selected cluster structures illustrating the core ar- 
rangements corresponding to various (cr,e) values are 
shown in Fig. O We see from Figs. ^ and that the 
X5Y2 cluster exhibits core X-atom structures that in- 
clude trigonal bipyramidal, planar and square pyramidal 
core geometries. Of these, only the trigonal bipyramidal 
form is stable in the parent X5 system. This illustrates 
that a suitable choice of the (cr,e) parameters can control- 
lably induce core geometries not present as stable minima 
in the bare cluster. For example, the square pyramid core 
structure, seen in Fig. 02 as a stable system, corresponds 
to a transition state in the bare X5 cluster. 

Fig. El represents the X5Y2 cluster at four points in 
Fig.Hdefined by the (a,e) coordinates (0.4,0.5), (0.4,1.0), 
(0.4,1.5) and (0.4,2.0). Here the pairs of coordinates cor- 
respond to (a), (b), (c) and (d) part of Fig. ^ respec- 
tively. In other words, we keep value of cr=0.4 fixed, 
while increasing the value of e. Each disconnect ivity 
graph shows all inherent structures available to the sys- 
tem for the given (cr,e) values. The global minimum of 
each system is labeled by number 1 and contains as a 
recognizable component the square pyramid core struc- 
ture (see Fig. ISI2). In Fig. El a the square pyramid core 
structure is connected to two inherent structures, labeled 
by 2 and 3, by pathways whose energies do not exceed 
— 13.8 (in units of e^^). Since isomer # 2 contains the 
same core structure (the square pyramid) as the global 
minimum the corresponding isomerization thus does not 
lead to a change in the core structure of the cluster. 



For present purposes, therefore, the barrier that con- 
nects them is not a "relevant" barrier. The relevant 
barriers are those that connect inherent structures that 
contain different core structures. The inherent structure 
# 3 contains as the core structure a (distorted) trigo- 
nal bipyramid (see Fig. Ol). Therefore, the isomeriza- 
tion barrier that connects the inherent structure # 3 
with global minimum is the lowest relevant isomeriza- 
tion barrier and its value is AEi^3=0.986e^^. Fig. (HI b 
and Fig.Elc show that increasing the value of e increases 
isomerization barriers that connect inherent structure #1 
(the square pyramid core structure) with inherent struc- 
ture #2 (the distorted trigonal bipyramid core struc- 
ture). Numerically, these barriers are AEi^2=l-227e^^ 
and AEi^2=l-431e^^, respectively. In Fig.jSld the square 
pyramid core structure is connected by an isomerization 
barrier of AEi ^2=1 -6476^^ with two (almost degenerate 
in energy) distorted trigonal bipyramid core structures. 
As illustrated in Fig.jSland discussed above, the barriers 
that determine the isomerization kinetics of these newly 
induced structures are sensitive to the (cr,e) values and 
can thus be at least partially controlled. These two sim- 
ple results are specific demonstrations of goals 1 and 3 
stated above. 




(5.4) (5.5) 




FIG. 5: Plots of X5Y2 structures for selected (cr, e) values. 
The decimal number for each figure denotes the corresponding 
(cr, e) domain in Fig. 0] 

As a second illustration, we consider mixed clusters 
of the type X7Y3. This system builds upon a parent. 
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FIG. 6: Disconnect ivity graph for X5Y2 (cr, e) values demonstrating that we can control barriers for the selected inherent 
structures. The energy scale is in units of . The (cr, e) values for panels (a-d) are (0.4,0.5), (0.4,1.0), (0.4,1.5) and (0.4,2.0), 
respectively. 




FIG. 7: Ecore(cr, e) for X7Y3. Format for the plot is the same 
as m Fig.H 



seven-atom, "magic number" system known to exhibit 
a set of four, energetically distinct inherent structures. 
The core inherent structures and associated energies for 
the stable X7 inherent structures are presented in Fig. [3 
Figure [3 a (a, e) contour plot of the core-atom poten- 
tial energies of the lowest total energy X7Y3 clusters, 
again reveals the presence of definite "core-phases" . As 
illustrated in Fig. |S1 some of these regions correspond to 
various core structures present in the parent X7 system 
while others correspond to new structures not seen in 
the original, single-component cluster. We can see from 
Figs. [7|-|ni that the impurity Y-atoms provide us with 
significant control over the relative ordering of the core 
energies of the parent X7 system. Specifically, by choos- 
ing an appropriate range of (a, e) values, we can generate 
X7Y3 clusters in which the lowest (total) energy inherent 
structure can have core structures that are either pentag- 
onal bipyramid, capped octahedral, or bicapped trigonal 
bipyramidal in nature. Moreover, since we can manipu- 
late the isomerization barriers in these systems, we can 
at least partially stabilize clusters that exhibit selected 
core structures with respect to isomerization. This is il- 
lustrated in Fig. [HI 

Figures El a - Eld represent the X7Y3 cluster at four 
points in Fig. [3 with X7Y3 (cr,e) coordinates (0.4,0.5), 
(0.4,1.0), (0.4,1.5) and (0.4,2.0), respectively. The num- 
ber of inherent structures available to the X7Y3 cluster 
varies from more than 800 in Fig.[Ula to 400 in Fig. Eld. 
Since we are primarily interested in energetically low- 
lying inherent structures we show only lowest 70 inher- 
ent structures. The global minimum of each system is 
labeled by number 1 and contains as a recognizable com- 
ponent the core structure shown in Fig. IHl5. We should 
mention that for a given range of a and e values, (a, e)G 
[0.1,2.0], we have not been able to find a global minimum 
that would contain as a recognizable component inherent 
structure # 4 of the parent X7 cluster (see Fig.|21d). This 




FIG. 8: Plots of selected X7Y3 structures for various (a, e) 
values identified in Fig. Q The number of the structures cor- 
respond to the regions labeled in Fig. |3 Note that many of 
the core structures for these systems are not stable energy 
structures of the bare X7 system. 



is the reason why none of the domains in Fig.Qis labeled 
by number 4. In Fig. Ela the global minimum, the core 
structure # 5 (see Fig.|Hl5), is linked to inherent structure 
# 2 which contains the (distorted) capped octahedron 
core structure (see Fig. |Hl2). The isomerization barrier 
between them is AEi,2=0.494e^^. Fig.Elb and Fig.Elc 
show that increasing the value of e increases isomerization 
barriers, that connect inherent structure #1 (see core 
structure in Fig|Hl5) with inherent structure #3 and #2 
(a distorted capped octahedron core structure), respec- 
tively. Numerically, these barriers are AEi^3=0.975e^^ 
and AEi^2=l-136e^^, respectively. The inherent struc- 
ture #2 in Fig. Elb contains the same core structure as 
the global minimum and, therefore, has not been con- 
sidered relevant for the isomerization (see explanation 
above for Fig. Ela). In Fig. Eld the core structure # 5 
is connected to two, energetically almost degenerate, in- 
herent structures labeled by 2 and 3, by pathways whose 
energies do not exceed —40.5 (in units of e^^). Simi- 
larly to the case of Fig. Ela, the isomer # 2 contains the 
same core structure as the global minimum so the cor- 




FIG. 9: Disconnectivity graph for X7Y3(cr, e) values demonstrating that we can control barriers for the selected inherent 
structures. The energy scale is in units of e^^ . Only branches leading to the 70 lowest-energy minima are shown. 
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responding isomerization does not lead to a change in 
the core structure of the cluster. The barrier that links 
them is not a relevant barrier. The inherent structure 
# 3 contains as the core structure a distorted capped 
octahedron. Therefore, the isomerization barrier which 
connects the inherent structure # 3 with global minimum 
is the lowest relevant isomerization barrier and its value 
is AEi,3=1.304e^^. 

IV. CONCLUSIONS 

In the present work, we have considered the general 
task of altering core cluster structures. We are, in effect, 
attempting to turn the logic of the minimization prob- 
lem upside down. Rather than seeking the global mini- 
mum of complex potential energy surfaces, we are instead 
attempting to exploit what has been learned about the 
general minimization problem to controUably alter core 
cluster structures. Specifically, we are examining the ex- 
tent to which we can induce new core geometries as well 
as reorder and stabilize existing, higher-lying, local core 
structures. 

Our approach, in the present discussion, has been ther- 
modynamic in nature. We have utilized selected adatoms 
to effect our desired core cluster modifications. We have 
presented results for two simple binary cluster examples, 
the X5Y2 and X7Y3 systems, to validate our approach. 



We speculate that there are at least two important di- 
rections for future theoretical development of the present 
ideas. One direction will be to explore the use of more 
complex adsorbates to achieve selected core cluster struc- 
tures. One could, for example, imagine using "exterior" 
methods in which encapsulating agents of well-defined 
geometries were utilized to induce desired core struc- 
tures. Alternatively, "interior" approaches in which com- 
plex objects, perhaps even previously engineered clusters, 
could be utilized as "seeds" or "templates" to produce a 
desired structure in the surrounding cluster (either glob- 
ally or locally). Another important direction will be to 
explore the extent to which previously engineered clus- 
ter structures can be assembled using "cluster assembled 
materials" methods to produce larger scale, macroscopic 
structures. If this proves possible, it would seem to offer 
an important direction in the production of novel mate- 
rials starting from synthetic precursors whose core struc- 
tures and properties are highly varied and are under user 
control. 
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